The characteristics of thermahzation of boost-invariant plasma from holography 
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We report on the approach towards the hydrodynamic regime of boost-invariant A/" = 4 super 
Yang-Mihs plasma at strong coupUng starting from various far-from-equihbrium states at r = 0. 
The results are obtained through numerical solution of Einstein's equations for the dual geometries, 
as described in detail in the companion article arXiv: 1203.0755 Despite the very rich far-from- 
equilibrium evolution, we find surprising regularities in the form of clear correlations between initial 
entropy and total produced entropy, as well as between initial entropy and the temperature at 
thermahzation, understood as the transition to a hydrodynamic description. For 29 different initial 
conditions that we consider, hydrodynamics turns out to be definitely applicable for proper times 
larger than 0.7 in units of inverse temperature at thermahzation. We observe a sizable anisotropy in 
the energy-momentum tensor at thermahzation, which is nevertheless entirely due to hydrodynamic 
effects. This suggests that effective thermahzation in heavy ion collisions may occur significantly 
earlier than true thermahzation. 



Introduction. One of the outstanding problems of the 
dynamics of quark-gluon plasma (QGP) is the under- 
standing of the physics of thermahzation. In relativistic 
heavy-ion collisions at RHIC and LHC the quantitative 
description of experimental data requires the applicabil- 
ity of hydrodynamics from a very early stage [1,. How- 
ever, our insight into the non-equilibrium dynamics of 
QGP is very scarce. The above problem is often referred 
to as 'the early thermahzation puzzle'. This is in fact a 
misnomer as viscous hydrodynamics may turn out to be 
applicable when the pressures are still quite anisotropic, 
going against the commonly accepted paradigm that true 
thermahzation is necessary. One of the main results of 
the present work is that for a wide range of initial con- 
ditions this is indeed the case. Subsequent isotropization 
towards true thermodynamic equilibrium occurs purely 
within the quantitatively well understood viscous hydro- 
dynamics and is trivial in comparison. 

The key physical question of interest is the time scale 
after which viscous hydrodynamic description becomes 
valid. This has a further refinement as viscous hydro- 
dynamics is really a gradient expansion with new trans- 
port coefficients appearing at each order. So it is very 
interesting to determine to what extent would all- order 
resummed hydrodynamics describe the plasma evolution 
and to what extent is one forced to incorporate genuine 
non- hydrodynamic degrees of freedom. Furthermore, the 
dynamics of plasma expansion will strongly depend on 
the initial state. It is very important to understand if 
there exists some simple physical characterization of the 
initial state determining the characteristics of the tran- 
sition to hydrodynamics and subsequent evolution. Fi- 
nally, it is interesting to understand the amount of en- 
tropy produced during different stages of the dynamics. 

In this letter we will address the above questions 



for plasma configurations invariant under longitudinal 
boosts and with no dependence on transverse coordi- 
nates. This kinematical regime was first introduced by 
Bjorken [2] and roughly mimicks an infinite energy colli- 
sion of infinitely large nuclei. 

Within QCD there are no techniques allowing to ad- 
dress these issues from first principles. It is thus quite 
natural to consider the same questions in the context 
of strongly coupled plasma in the A/" = 4 supersymmet- 
ric gauge theory for the description of which one can 
use the AdS/CFT correspondence [3 . There, the time- 
dependence of plasma is translated into gravitational dy- 
namics in 5 dimensions with a negative cosmological con- 
stant and appropriate boundary conditions. Using these 
methods perfect fluid hydrodynamics was derived at the 
nonlinear level in the boost-invariant setting [4J , the value 
of shear viscosity was shown to agree [5 with the one ex- 
tracted from linear perturbations [6], and finally viscous 
hydrodynamics was derived without any symmetry as- 
sumptions [7 . 

Once we consider the far-from-equilibrium regime for 
small proper times, gradient or scaling expansions cease 
to be valid, and one has to deal with full Einstein's 
equations. Previous work by some of us [8 , motivated 
by the early results of [9 , used power series expansions 
around r = to study strongly non-equilibrium regime 
of Bjorken flow. Unfortunately, the radius of convergence 
of these power series was insufficient to analyze the tran- 
sition to hydrodynamics. On the other hand, the numer- 
ical work of [10 necessarily introduced a deformation of 
the physical 4-dimensional metric to pump energy and 
momentum into the vacuum at early times and create 
in this way a far-from-equilibrium state. Such a way of 
generating the initial state precludes the analysis of the 
physical evolution starting from r = 0, in particular the 
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investigation of the influence of the initial conditions on 
thermahzation and entropy production that we are inter- 
ested in. 

Motivated by this, we developed a new numerical 
framework using the ADM formalism of numerical rel- 
ativity and analyzed the evolution of the plasma sys- 
tem starting from a range of initial conditions. These 
correspond, in our setup, to specifying a single metric 
coeflicient function ( 'initial profile ^) for the initial geom- 
etry on the hypersurface r = 0. The initial hypersur- 
face is the same as in [8 , however without any spuri- 
ous coordinate singularities. Subsequently we solve nu- 
merically 5-dimensional Einstein's equations and obtain 
plasma energy-momentum tensor from the asymptotics 
of the solution at the AdS boundary. The details of this 
setup can be found in a companion article pjj, while 
in the present letter we will concentrate on the physical 
questions mentioned above. 

Boost- invariant plasma and hydrodynamics. The 
traceless and conserved energy-momentum tensor of a 
boost-invariant conformal plasma system with no trans- 
verse coordinate dependence is uniquely determined in 
terms of a single function (Trr) - the energy density at 
mid-rapidity e{r). The longitudinal and transverse pres- 
sure are consequently given by 

d ^ Id,. 

PL = -£-r--e and pT=e^-T--e. (1) 
dr 2 dr 

It is quite convenient to eliminate explicit dependence 
on the number of colors Nc and degrees of freedom by 
introducing an effective temperature Tg// through 
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Let us emphasize that Tg// does not imply in any way 
thermahzation. It just measures the temperature of a 
thermal system with an identical energy density as €{r). 

All order viscous hydrodynamics amounts to present- 
ing the energy-momentum tensor as a series of terms ex- 
pressed in terms of flow velocities and their deriva- 
tives with coeflicients being proportional to appropriate 
powers of Tg//, the proportionality constants being the 
transport coefficients. For the case of A/" = 4 plasma, 
the above mentioned form of T^j^ is not an assumption 
but a result of a derivation from AdS/CFT [7]. Hydro- 
dynamic equations are just the conservation equations 
d^T^^ = 0, which are by construction first-order differ- 
ential equations for Tg//. 

In the case of boost-invariant conformal plasma this 
leads to a universal form of first order dynamical equa- 
tions for the scale invariant quantity w = T^f f ■ r namely 



r d 
w dr 
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hydro 



(w) 
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(3) 



where Fhydro{u)) is completely determined in terms of the 
transport coefficients of the theory, much in the spirit of 
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FIG. 1. a) F(w)/w versus w for all 29 initial data, b) Pressure 
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for a selected profile. Red, blue and green 



curves represent 1**, 2^ and 3^ order hydrodynamics fit. 

\TT . For A/" = 4 plasma at strong coupling Fhydro{uj) / w 
is known explicitly up to terms corresponding to 3^^ order 
hydrodynamics |^ 



1 - log 2 15 - 27r2 - 45 log 2 + 24 log^ 2 
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(4) 

The importance of formula ([3| lies in the fact that if the 
plasma dynamics would be governed entirely by (even 
resummed) hydrodynamics including dissipative terms 
of arbitrarily high degree, then on a plot of ^£fW = 
F{w)/w as a function of w trajectories for all initial con- 
ditions would lie on a single curve given by Fhydroi'^)/'^- 
If, on the other hand, genuine non-equilibrium processes 
would intervene we would observe a wide range of curves 
which would merge for sufficiently large w when thermal- 
ization and transition to hydrodynamics would occur. 

In Figure la we present this plot for 29 trajectories cor- 
responding to different initial states. It is clear from the 
plot that non-hydrodynamic modes are very important in 
the initial stage of plasma evolution, yet for all the sets 
of initial data, for w > 0.7 the curves merge into a single 
curve characteristic of hydrodynamics. In Figure lb we 
show a plot of pressure anisotropy 1 — = 12^^^^ — 8 
for a selected profile and compare this with the corre- 
sponding curves for 1^^, 2^^ and 3^^ order hydrodynam- 
ics. We observe, on the one hand, a perfect agreement 
with hydrodynamics for w > 0.63 and, on the other hand, 
a quite sizable pressure anisotropy in that regime which 
is nevertheless completely explained by dissipative hy- 
drodynamics (see [10^ for similar conclusion). 

In order to study the transition to hydrodynamics in 
more detail, we will adopt a numerical criterion for ther- 
mahzation which is the deviation of r:^^; from the 3^^ 
order hydro expression ^ 
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< 0.005. 
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Despite the bewildering variety of the non-equilibrium 
evolution, we will show below that there exist, however, 
some surprising regularities in the dynamics. 
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FIG. 2. The apparent horizon (black u-shaped curve) and a 
radial null geodesic (red curve) sent from the boundary (left 
edge of the plot) at r = into the bulk for a sample profile. 



Initial and final entropy. Apart from the energy- 
momentum tensor components, a very important physi- 
cal property of the evolving plasma system is its entropy 
density s (per transverse area and unit (spacetime) ra- 
pidity). In the general time-dependent case, the precise 
holographic dictionary for determining entropy is miss- 
ing. Nevertheless in the present case due to high sym- 
metry, entropy seems to be defined unambiguously in 
terms of I/AGn of the apparent horizon (AH) area el- 
ement mapped onto the boundary along ingoing radial 
null geodesies [IQl IT4l 115] . 

For all of the initial profiles that we considered we ob- 
served an apparent horizon which was pierced by a radial 
null geodesic starting from r = on the boundary (see 
Figure 2). This shows that the initial conditions had 
always some entropy per unit rapidity to start with. 

The main very surprising observation of our work is 
that the initial entropy density measured in units of ef- 
fective temperature at r = is a key characterization of 
the initial state which, to a large extent, determines the 
features of the subsequent transition to hydrodynamics 
as well as the final produced entropy. 

In the following it is convenient to introduce a dimen- 

sionless entropy density Sn-eq = sah/ (^iN^Tr'^iT^'ff)'^^ 
In order to evaluate the final entropy density at r = oo, 
we adopted the following strategy. After observing a pas- 
sage to hydrodynamics, we fitted 3^^^ order hydrodynamic 
expression for Tg// 
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to obtain the remaining single scale A. Since at r = oo 
perfect fluid hydrodynamics applies, we can use the stan- 
dard expression for entropy to get s^fl^^ = A^ • (T^^ff) 

Once this has been done we can now determine the 
entropy production s^^l^^ — s^^l^^ as a function of 
for all the considered profiles. Despite the huge differ- 
ences in the evolution evident in Figure la, we observe 
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FIG. 3. Entropy production as a function of initial entropy. 



a clear functional dependence of the entropy production 
on the initial entropy. The results are shown in Figure 3 
together with a fit s\{\q - = 1-64 • {s^n-eqf'^^' 

Properties of thermalization. We will now proceed 
to study in detail the properties of the transition from 
far-from-equilibrium regime to hydrodynamics. We will 
adopt the criterion ([5|, which imposes quite precise 
agreement between the equations of motion coming from 
third order hydrodynamics (being the most precise de- 
scription currently known) and the actual evolution of 
the energy density of the plasma obtained from numeri- 
cally solving the full Einstein's equations. This criterion 
is quite different from criterions based on isotropization 
of the longitudinal and transverse pressures like the one 
adopted in [8 . In particular. Figure lb shows quite a 
sizable pressure anisotropy, which is nevertheless entirely 
due to hydrodynamic modes. 

Using the criterion ([5|, we determine the thermaliza- 
tion times for 29 initial profiles. If we were to modify the 
threshold, the thermalization time would of course shift 
but for most profiles not significantly ^J. However, it is 
fair to say that thermalization is not a clear-cut event but 
rather happens in some narrow range of proper times. 

With this proviso we will now proceed to analyze the 
following features of the thermalization time: (i) the di- 
mensionless parameter w = Tg// • r, (ii) the thermaliza- 
tion time in units of initial temperature and (iii) the ratio 
of the effective temperature at the time of thermalization 
to the initial (effective) temperature. 

In Figure 4, we show a plot of the values of w at the 
time of thermalization as a function of the initial en- 
tropy. We see that for a wide range of initial entropies, 
the values of w at thermalization are approximately con- 
stant and decrease only for initial data with very small 
entropies. 

Subsequently, we found an unexpectedly strong corre- 
lation of the thermalization time with the initial entropy 
(see Figure 5). This is very surprising taking into ac- 
count the huge qualitative differences in the evolution of 
the plasma when starting from various initial conditions. 

Another important aspect is the question which part 
of the cooling process of the plasma occurs in the far- 
from-equilibrium regime and which part occurs within 
hydrodynamic evolution. This can be quantified by the 
ratio of the effective temperatures at thermalization time 
and at r = 0. The plot in Figure 5 shows very clear 
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FIG. 4. The dimensionless parameter w — rTeff at thermal- 
ization as a function of the initial entropy s^^lg^. The straight 
hne corresponds to w^*^^ = 0.7. 
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FIG. 5. The thermaUzation time in the units of initial effective 
temperature (left) and the ratio of thermaUzation and initial 
effective temperatures T^^^^ jT^^^^ (right) as functions of the 
initial entropy . 



An effective thermalization time w^^^^ = T^^^V^^^^ i.e. 
thermalization time measured in units of effective tem- 
perature at thermalization depends on the initial state, 
but not strongly, and is between 0.37 and 0.67 for all 
considered initial states (a reasonable RHIC estimate 
T 500 MeV and r 0.25 /m/c gives w 0.63). Fi- 
nally, let us note that we could associate with all these 
initial data, an initial entropy already at r = due to 
the presence of an apparent horizon. This observation 
shows that thermalization and applicability of (all-order) 
viscous hydrodynamics is not necessarily associated with 
the sudden appearance of a horizon. 
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correlation of this quantity with the initial entropy. The 
meaning of the points with high entropy requires some 
comment. We found that for these initial conditions, 
the energy density initially rises and only later decreases, 
thus even a ratio of T^^ff^ / ^i// close to 1 is realized after 
a sizable non-equilibrium evolution pT]. 

Conclusions. The crucial new feature of the holo- 
graphic studies of Bjorken flow reported here is the 
ability to track physical observables from the far-from- 
equilibrium regime at r = up to thermalization and 
subsequent hydrodynamic evolution without introducing 
any deformations in the field theory lagrangian. The ini- 
tial state is highly anisotropic, in particular always has 
a negative longitudinal pressure [8, 9 . Despite the very 
rich early time dynamics, which depending on the ini- 
tial state might have a plateau, a bump or a sharp de- 
crease in the effective temperature as a function of proper 
time, we uncovered surprising regularities in the behav- 
ior of total produced entropy and effective temperature 
at thermalization as functions of initial entropy (all mea- 
sured in units of effective temperature at r = 0). An 
interesting curiosity is that despite describing an expand- 
ing medium, the effective temperature at thermalization 
might be higher than the initial one for initial states 
with sufficiently big entropy. For initial states with small 
entropy, the energy density at thermalization is much 
smaller than the one at r = 0, and hence a significant 
part of the cooling process is of a non-equillibrium na- 
ture. Moreover we observe generically a sizable pressure 
anisotropy at thermalization, which is nevertheless en- 
tirely understood in terms of dissipative hydrodynamics. 
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